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ABSTRACT 

This paper introduces the formahsm which connects between rotation measure (RM) 
measurements for extragalactic sources and the cosmological magnetic field power spec- 
trum. It is shown that the amphtude and shape of the cosmological magnetic field power 
spectrum can be constrained by using a few hundred radio sources, for which Faraday 
RMs are available. This constraint is of the form B^ms 1 x [2.6 x 10^^cm~'^/nf,]/i 
nano-Gauss (nC) on ~ 10 — 50 h~^Mpc scales, with fii, the average baryon density and 
h the Hubble parameter in 100 km s~^ Mpc~^ units. The constraint is superior to and 
supersedes any other constraint which come from either CMB fluctuations, Baryonic 
nucleosynthesis, or the first two multipoles of the magnetic field expansion. The most 
adequate method for the constraint calculation uses the Bayesian approach to the max- 
imum likelihood function. I demonstrate the ability to detect such magnetic fields by 
constructing simulations of the field and mimicking observations. This procedure also 
provides error estimates for the derived quantities. 

The two main noise contributions due to the Galactic RM and the internal RM 
are treated in a statistical way following an evaluation of their distribution. For a 
range of magnetic field power spectra with power indices —1 < n < lina flat 
cosmology (17^=1) we estimate the signal-to-noise ratio, Q, for limits on the mag- 
netic field Brms on ~ 50 h^^Mpc scale. Employing one patch of a few square de- 
grees on the sky with source number density Usrc, an approximate estimate yields 
Q ~ 3 X {Brms/lnG){nsrc/50deg~'^){2.6 x lO^'^cm-V^b) ^- An all sky coverage, with 
much sparser, but carefully tailored sample of ~ 500 sources, yields Q ~ 1 with the 
same scaling. An ideal combination of small densely sampled patches and sparse all-sky 
coverage yields Q — 3 with better constraints for the power index. All of these estimates 
are corroborated by the simulations. 

Subject headings: Cosmology: Large-Scale Structure of Universe, Theory - Magnetic 
Fields - Polarization - Methods: Statistical 
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1. 



INTRODUCTION 



There is plenty evidence for magnetic fields on the scale of the earth, the solar system, the 
interstellar medium, galaxies and clusters of galaxies [for the extragalactic magnetic field see review 
by Kronberg (1994, K94) and references therein]. We still don't know however whether any larger 
scale magnetic fields exist. There are various upper limits on the rms value of the cosmological 
magnetic field, but the upper limits are still too high to exclude any relevance to magnetic fields 
observed on smaller scales. Cosmological magnetic fields might help to resolve the bothersome 
question of the origin of the ~ 10~^G magnetic field which exists on galactic scales or fields of 
similar size seen on the scales of galaxy clusters. A primordial (pre - galaxy formation) magnetic 
field amplified by the gravitational collapse, and then by processes like differential rotation and 
dynamo amplification (Parker 1979), can serve to seed these observed smaller scale magnetic fields. 

The most stringent upper limits to date for magnetic fields on cosmological scales come from 
three different sources: 

1. Big Bang Baryonic Nucleosynthesis (BBN). 

Magnetic fields that existed during the BBN epoch would affect the expansion rate, the 
reaction rates, the electron density (and possibly the space-time geometry). By taking all 
these effects into account in the calculation of the element abundances, and then comparing 
the results with the observed abundances one can set limits on the magnetic field amplitude. 
The limits for homogeneous magnetic fields on scales 3> 10'^"^ h^^Mpc (distances are quoted 
in comoving coordinates) and up to the BBN horizon size (~ 10~^ h~^Mpc) are 
G in terms of today's values (Grasso & Rubinstein 1995, 1996 ; Cheng, Olinto, Schramm, & 
Truran 1996). The relevance of these limits for the maximum value of magnetic field seeds on 
the subgalactic scale is obvious, but in order for these limits to be relevant to the intergalactic 
magnetic field, further assumptions about the super-horizon magnetic field power spectrum 
and the magnetic field generation epoch should be made. 

2. The CMB radiation. 

The presence of magnetic fields during the time of decoupling causes anisotropy, i. e. different 
expansion rates in different directions (Zel'dovich &; Novikov 1975). Masden (1989) provides 
a useful formula for the connection between the measured temperature fluctuations on a given 
scale at z^,, the decoupling time, and the limits on the equivalent scale today. For a magnetic 
field energy density much smaller than the matter energy density, the limit is 



For typical values of /i = 1 (the Hubble constant in units of 100 km s"^ Mpc"^), temperature 
fluctuations of ~ IQ-^, J7 = 1 and — 10^ the limit becomes Brms < 10-^ — lO-^G. Recently, 
Barrow, Frreira, Sz Silk (1997) used the COBE 4-year data to constrain the magnetic field on 
the horizon scale by B{R = Rh) < 6.8{nh?)^/^ nG. 




(1) 
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3. Multipole expansion of RM observations. 

Passing through a magnetoionic medium, polarized light undergoes rotation of the polariza- 
tion vector (Faraday rotation, for definition see §^). The rotation amount is proportional to 
the dot product of the magnetic field and the light propagation direction, and to the distance 
the light travels through the medium. If any cosmological magnetic field exists, the farther 
away the source, the more the polarized light component is rotated. We refer to this as the 
"monopole" term. This term exists even if the field has no preferred direction on the horizon 
scale, namely even without breaking the isotropy hypothesis. 

If there is a preferred direction to the field on a cosmological scale, a dipole will show up in 
the RM measurements across the sky, due to the different sign of the dot product (and thus 
the different rotation direction). So far attempts have been made only to identify a monopole 
or a dipole in the magnetic field for scales of 2; ~ 2.5 (Kronberg 1976) and z ~ 3.6 (Vallee 
1990). The search for a dipole signature in the Faraday rotation values for a sample of extra- 
galactic sources (QSO) yields limits of i3 < 1 — 10^^ nG depending on the assumed cosmology 
(Kronberg & Simard-Normandin 1976, Kronberg 1976). Vallee (1990) tried to estimate the 
limits on the dipole and the monopole terms and concluded that an all-prevailing (up to 



3.6) field must be smaller than 6 x 10 



10" 



nG . 



Other potential methods of detecting a possible cosmological magnetic field include distortion of the 
acoustic ("Doppler") peak in the CMB power spectrum (Adams, Danielsson, Grasso, & Rubinstein 
1996) and Faraday rotation of the CMB radiation polarized components (Loeb & Kosowsky 1996). 
These methods where not designed to probe any magnetic field generated after recombination. 
Their implementation is still pending on the upcoming measurements of the CMB fluctuations and 
polarization. 

A method that does probe magnetic fields generated in the post recombination era, relies on 
cosmic ray (CR) detection. The effect of magnetic fields on high energy CRs (> 10^^ eV) is twofold. 
It will first alter the energy distribution of CRs (Lee, Olinto, & Sigl 1995, Waxman & Miralda- 
Escude 1996) and then, if the CR direction is identified and attributed to a known region or source, 
the magnetic field either smears the directionality (for field coherent scales much smaller than the 
distance to the source) or deflects the CR direction from aligning with the source direction. Current 
limits from this last effect are either very weak or non-existing. 

Theories for primordial magnetic fields in the framework of structure evolution paradigms 
have been suggested by a number of authors. Magnetic fields on cosmological scales emerge in the 
framework of these theories in one of three epochs: the infiation era, the plasma era, and the post- 
recombination era. In the inflation era, magnetic fields form due to quantum mechanical processes 
(Turner & Widrow 1988; Quashnock, Loeb, & Spergel 1989, Vachaspati 1991, Ratra 1992a, 1992b, 
Dolgov & Silk 1993) or possibly magneto-hydrodynamics (Brandenburg, Enquist, & Olesen 1996). 
During the plasma era, prior to recombination, magnetic fields form due to vorticity caused by 
the mass difference between the electron and the proton (Harrison 1973, but see Rees 1987), or 
where magnetic field fluctuation survival (no vorticity assumed) depends on the fluctuation scale 
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(Tajima et al. 1992). Other authors have pointed out that even in the post recombination era 
(but before galaxy formation) magnetic fields can be generated by tidal torques (Zweibel 1988), or 
by cosmic strings wakes (Ostriker & Thompson 1987, Thompson 1990). Magnetic field generated 
in the process of structure formation itself, due to falling matter (Pudritz &; silk 1989), or star 
burst regions, are less likely to feed back into the intergalactic space and to affect the cosmological 
magnetic field. Nevertheless, proposals in this spirit were also considered by taking into account 
wind driven plasma which fills intergalactic space with magnetic fields. 

Whatever the origin of cosmological magnetic fields, each suggested theory, specifies the mag- 
netic field amplitude and power spectrum shape. Thus, each of the theories can in principle become 
refutable if we could measure, or put limits on, the normalized power spectrum of the cosmological 
magnetic field. 

This challenge and the realization that with existing or upcoming data, it will be possible to 
measure the magnetic field power spectrum has led us to develop a method for doing so. We argue 
that much better estimates (or upper limits) for the cosmological magnetic field can be calculated 
by using the RM correlation matrix. 

All previous attempts to limit the cosmological magnetic field dipole on very large scales 
(z ~ 2.5 — 3.6), suffer from two drawbacks. From a theoretical point of view it is difficult to 
reconcile the necessary anisotropy such a dipole imposes with other evidence which suggest global 
isotropy. From a practical point of view, this test is confined to one scale, and doesn't allow 
measurements over a whole range. Only by extending the measurement to a whole range, can one 
estimate the full power spectrum (PS). Although contribution from a "random walk" can be useful 
even if only a monopole is calculated (see below), the information is partial, inferior to the full 
derivation of the magnetic field by RM correlations, and sensitive to evolutionary effects. 

Why is the correlation approach more advantageous than other type of statistics and in partic- 
ular why is it better than the monopole approach? In the monopole approach we are looking for a 
correlation between RM values and the source redshift. Let's pretend at first that all measurements 
are exact and there are no noise sources (we'll notice in §^ that these two assumptions are far too 
optimistic). Consider a scale over which the magnetic field is coherent with an rms value of 
Byms- A line-of-sight to a source at a distance rgrc will cross A'^; = r^rc/^o such regions. There is no 
correlation in the magnetic field orientation between the regions, and on average each one of them 
contributes RM*„^ oc B„nsh- The overall contribution due to A^; random walk steps sums up to 
RMrms = ^f^i^^^rms ^ BrmsVloi^src ■ It is thus clear that the larger rgrc, the bigger RM^ms we 
expect for an ensemble of sources located at rgrc ■ Now let's turn to the correlation approach. Two 
lines-of-sight to two adjacent sources (separation angle 7 and assume the same redshift for the two) 
are observed for correlation between the RM values of the polarized light emanating from them. 
If 'jrsrc ^ Iq the two light rays undergo the same rotation in each patch of length Iq (for small 7 
we neglect the difference in the cos{9) term between the magnetic field and the line-of-sight). The 
ensemble average of tli6 correlation term would, then be <C R,]V[iR,]Vl2 ^^^^ — A^/FtlVI^^^ cx BY-ms^src-^ 
with the same proportionality constant as before. In the limit of lo Vsrc, the two expressions 
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coincide, but for any value of Iq < Vgrc, the signal from the correlation calculation is amplified 
by the factor ^r^rc/^o- For example, if Vgrc = 1700 h^^Mpc {z ~ 1) and Iq = 10 h~^Mpc the 
signal is amplified by more than an order of magnitude. Moreover, the correlation provides infor- 
mation about the power spectrum (different 7 and r^rc values) of the magnetic field that otherwise 
we wouldn't be able to obtain. The correlation also makes it easyer to separate the noise from 
the signal, unlike the monopole approach. We will return to this approximation later on, under 
more realistic considerations (§3^), when we attempt to evaluate the number of pairs needed in 
order to establish a certain signal-to-noise ratio at a given scale (7). For now, this simple model 
demonstrates the advantages of this paper's approach. 

In order to relate actual data of RM measurements to the magnetic field, we begin by con- 
sidering the connection between the two. In section ^ we introduce the magnetic field correlation 
tensor and its connection to the RM correlation. In section |^ we make use of these definitions and 
describe the calculation procedure for the RM correlation from the raw data to the constraints on 
the magnetic field power spectrum. A key issue is the estimate of noise from non-cosmological con- 
tributions to the RM. The impatient reader can turn immediately to §^]^ to get a rough estimate 
for the expected signal-to-noise ratio. We demonstrate the procedure in the following section (§4), 
where we simulate a few realistic examples of cosmological magnetic fields, and exploit them to get 
RM measurements from which we derive back an estimate for the original power spectrum. The 
simulations also allow us to perform a realistic error analysis. In section 5 we discuss prospects for 
applying the procedure to real data and conclude with our results. 



2. THE ROTATION MEASURE CORRELATION MATRIX 

Linearly polarized electro-magnetic radiation of frequency v traveling a distance dl through 
a non-relativistic plasma medium with the magnetic field B, is rotated by the angle d(/>. In other 
words, the polarization angle is changed by (e.g. Lang 1978 Eqs. 1-268 — 1-270) 



27rm^z/^ 



B-dl, (2) 



where e, me, and Ug are the electron charge, mass, and number density respectively (we use units 
of c = 1). The rotated polarization vector itself is written 

Kz.) = |p(l.)|e2^^ (3) 

a pseudo- vector degenerated in rotation of (j) by nvr. It is related to the Stokes parameters via 

\P\ = J ; ^=2*-^y> (4) 

that in turn are expressed by the electric field components perpendicular to the wave propagation 
direction (z) 

I = El + El 



- 6 - 



Q - E^- Ey 

U = E;E., + EyEl (5) 
where time average is explicitly assumed. 

In the cosmological context, some of the quantities are functions of time (or equivalently redshift or 
distance). We assume fully ionized gas (ionization fraction X(, = 1) and thus the average number 
density of free electrons is 

fie = XeUh = fib ■ (6) 

The baryon number density at a specific location is ni){x) = n{y{l+5{x)), with 5{x) the dimensionless 
density fluctuation. Number densities scale like a~^, where a is the scale factor. The frequency 
scales like a~^, and under the assumption of flux conservation, B scales like a~^. The distance unit 
"dr' in a Robertson- Walker metric is given by 

dr 

dl = dt = a(t)d7? = _ , (7) 

expressed via r] the conformal time. Since a = (1 + z)~^ for the redshift z, the overall rotation of 
the polarized radiation of a source at the redshift z, the direction q as observed from x, (we set 
X = for simplicity) and at frequency z/q, is given by 

1^ / [l + 5{qr')] {l+z'fBo{qr')-q = 8 

Today's values are all denoted by the subscript "0" and z' = z{r'). The time dependence of 5{qr') 
is taken into account later on (c/. Eq. 12). The two point correlation function of the RM is defined 
by 

T = {RMiqi,zi)RM{q2,Z2)), (9) 

see Nissen & Thielheim (1975) for a similar definition. Throughout the paper, (...) is the notation 
for an ensemble average. We use the notation f ' = ri ' — r2 ' and obtain 

2 



/ , ' , [l + 5{q,r[)]{l + z[f ^ [1 + 5{q,r',)] (1 + 



J ^3 



X 9ii^2jY / B^i{x)Boj{x^r')i?x, (10) 

where summation over the spatial components i,j is assumed. The integral (^) is made out 
of four terms. Two of the terms (those involve only one power of S) vanish because ((5) = 
and we assume vanishing correlation between the density fluctuations and magnetic fluctuations 
\{B{x)b{x + r)) = OV(x, r)] due to the vector nature of B and the scalar b. 

The two remaining terms are 

Jo Ji _ Kr? Jo . i _ Kr'^ 
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QiiQ2j^ J Boi{x)Boj{x + f')d^x, (11) 

and 

\27Tm^J Jo ^l^K^ - Kr'i 

qiihj^ J Bo,{x)Boj{x + r')d^x, (12) 

where ^(r) is the baryonic matter correlation function that we identify with the matter correlation 
function. Today's value of it is and the missing powers of (1 + z) in Eq. ( p!^ ) account for the 
linear evolution of this correlation. The RM correlation function is the sum T = Tq + T^. 

The last integral of Eqs. (IC, 11, ^) is the correlation tensor of the magnetic field defined as 
Cij{r) = {Bi{x)Bj{x + r)). In the cosmological case, we assume isotropy on scales much smaller 
than the horizon and confine Cij to be a function of |r| only. Being a divergence- free field the 
magnetic field correlation can be written (Monin & Yaglom 1975) as a combination of parallel (Cy ) 
and perpendicular (Cj_) functions. The parallelism and orthogonality are given with respect to the 
connecting vector: f = fi — ^2 

C,i{k) = [q|(fc) - C^{k)]^ + C^{k)6fi , (13) 

where S^j is the Kronecker (^-function. In the spectral domain we define the function E{k) = 
4:TTk'^C_i_{k) and express both correlation functions by the function E(k), namely 

C±{r) = j^dkE{k) [jo{kr) - ^-^^ ; C\\{r) = 2 dkE{k)^-^ , (14) 

with ji denoting the spherical Bessel function of the i^^ order. When people refer to the "magnetic 
power spectrum" they usually mean C^{k) which should further be multiplied by k'^ in order to get 
E{k). We shall hereafter comply with this notation and refer to C±{k) as the three-dimensional 
"magnetic power spectrum", Ps^k). Figure |^ show the two functions C||(r) and C±{r) in units 
of C±{0) for three values of the power index n = —1,0, 1 [PB{k) oc /c"]. Each function is shown 
with a minimal Gaussian smoothing of 1.5 h^^Mpc. Naturally, smaller power indices mean larger 
correlation length. 

Going back to the expression (|To| ) for the RM correlation we can rewrite it as 

^l-Kr? Jo Jl-Kr'f 



\2TTmlj Jo y^i _ Kr[^ 

{[qi(r') - Cx(r')](gi • r')(g2 • r') + Cx(r')gi ■ h} ■ (15) 



-8- 




Fig. 1. — The parallel (dotted-dashed line) and perpendicular (solid line) components of the magnetic field 
correlation tensor. Units are of Cj_(= Cy) in the origin. Plotted are three cases: (a) Pb(A;) oc k~^, (b) 
-Pb(^) = const., (c) Psik) oc k. Gaussian smoothing of 1.5 h~^Mpc applied. 

As expected, the RM correlation depends only on the relative distance, |r|, and the angle. 

For comoving distances between two sources at zi and Z2 separated by the angle 7 on the sky 
we use Weinberg's (1972, Eq. 14.2.7) coordinate transformation, implemented by Osmer (1981) 
for the qo = case (go - the deceleration parameter) and generalized by Boyle (1986) for other qo 
values. For the flat universe case (go = 0.5) the comoving relative distance reduces to the simple 
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form 



1 



1/2 



r{zi,Z2,j) = — ri{zi) +r2{z2) -2ri{zi)r2{z2)cosj . (16) 
-no L J 

The distance to a source at the redshift z (no cosmological constant) is given by (Kolb & Turner 

1990, Eq. 3.112) 

^ ^_ 1 21^oz + (21^o-4)(V17oz + l-l) 

"^''~Ho nlii + z) ' 

and can be easily generalized for other cases. 

In order to estimate the importance of [Eq. (p!2|)] we adopt the APM power spectrum 
(Baugh & Efstathiou 1993; Tadros & Efstathiou 1995) with the fitting formula 

P (A;) - ^ (^/^Q)'"'" (18) 

the fitting parameters m = 1.4, kc = 0.02, kQ = 0.19, and Gaussian smoothing of 1.5 h~^Mpc. 
We then derive the real space correlation function by the Fourier transform of the smoothed power 
spectrum. In the range 3 — 30 h~^Mpc the result is very similar to the derived APM real space 
correlation function (no explicit smoothing) = (?'/5.25 h~^Mpc)~^'^ (Baugh 1995). For this 
choice we get (regctrdless of the -Pb(^) forni, and. for Zg'p^ ^ 0.5), ITj/T] < 10^1 We hence identify 
T = To hereafter. 



2.1. The Power Spectrum Normalization 

We relate the magnetic field normalization to the power spectrum by Brms{R) where 

iB^J/)\R) = —j^ PB{k)k^WlTHikR)dk, (19) 

the factor 3 is due to the definition of Psik) [as C±{r)], and WG,TH{kR) is the Gaussian (Top-hat) 
window function of radius R, in k space. 



Wg = exp 



(kRY' 



Wth = jT^[sin{kR) - kRcos{kR)] . (20) 



A good benchmark to use for the normalization range is the CMB limit due to its simplicity and 
the availability of AT/T measurements on many scales (either today or in the near future). We 
therefore work in the range of ~ 1 nG throughout this paper. 

The observational limits on the magnetic field as deduced from the RM dipole [R ~ r(z = 2.5) 
or i? ~ r(z = 3.6)] are thus very different from the limits imposed by CMB fluctuations on 1' scale. 



We note that by Eq. (19) one cannot infer limits on the magnetic field magnitude from one scale 
to another. The limits must involve the power spectrum shape. Since no estimate for the latter 
exists, and since theoretical predictions for it range from power index ofn = — 3ton = 2 (Ratra 
1992a), limits on one particular scale hardly limit other scales. 
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Fig. 2. — RM correlation as function of the separation angle between source pairs at three different redshifts. 
The redshift values are Zg^.^ = 0.5, 1, 2 with thicker lines for higher redshift. The magnetic field normalization 
in all cases is 5^^(50 h~^Mpc) = 1 nG . Three magnetic power spectra are shown with the power indices 
n = —1, 0, 1. Note the different y axis scales. 

Figure |^ show the values for the RM correlation function for pairs of sources at the same 
redshift, and three different spectral indices. The normalization of the magnetic field in each case 
is i?^^(50 h~^Mpc) = 1 nG (Sl^ = 1, /i = 1, 0^ = 0.024). This scale resembles the mean separation 
between rich clusters. 

We notice that for close pairs, high redshift sources are preferable in order to get high signal, 
while a transition typically occurs at a fraction of a degree where it becomes preferable to use lower 
redshift sources. This transition is due to the uncorrelated magnetic fields experienced by high 
redshift sources with large separation. Only when the two lines-of-sight approach the correlation 
length, does contribution to the RM correlation becomes substantial. The different, uncorrelated 
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magnetic fields up to the correlation length, play the role of uncorrelated noise, and dominate the 
correlated signal from the low redshift segment of the line-of-sight. 

On the scale 0.5° — 2.5°, the RM correlation signal still stands on ~ 10— 10^ rad^ m~^ depending 
on the power spectrum. 



3. CALCULATION PROCEDURE FOR THE CORRELATION MATRIX 

3.1. The Raw Data 

The raw data consist of Ngrc extra Galactic objects for which polarization measurements are 
available in Nx wavelengths. Broten, Macleod, & Vallee (1988), and Oren &: Wolfe (1995) have 
emphasized the need for a careful selection of wavelengths for each observed source. For the 
current application there are additional special requirements that will become clear when we get 
to the discussion (§|5|). 

Each of such N\ measurements consists of two types of data : the degree of polarization, and the 
polarization position angle and 2(j) of Eq. (^]. Two functions are then constructed which give 
the dependence of these on wavelength. The RM is computed from the latter by using Eq. (^) where 
RM = (pX"'^. Each measured RM^ value has an error, 6^7^, that emerges from the measurement 
error (instrumental, ionosphere model, the earth magnetic field model etc.) and translates to the 
error in the fit from which the RM is derived. For typical sources with polarization degree of 
~ 10% (namely 10~^ — 10~^ Jy of polarized component at a few GHz frequency), the measurement 
error in each datum is of the order of a few degrees {e.g. Kato, Tabara, Inoue & Aizu 1987; 
Simard-Normandin, Kronberg, Button, 1980 (SKB); Oren &: Wolfe 1995). For realistic Nx = 4, 
this instrumental error translates to a typical error in the (f) — )? fit of 0.5 — 5radm~^ and is dwarfed 



by other noise terms in the procedure (c/. §3.2, §^^). We therefore set — from now on. 

The measured RMs are the sum of a few contributions, and cannot be directly plugged into 
the expression (^) to evaluate T. 

To begin with, there exists the internal RM/ for every one of the A^src extragalactic sources. Then 
there is the integrated RMc which is presumably due to the cosmological magnetic field (and 
the free electrons) along the line-of-sight to the source. There may also exist contribution RMj 
from intervening systems along the line-of-sight (foreground screen) either next to the source itself 
{e.g. the host galaxy of a quasar), Lyman-a systems, galaxies, or clusters of galaxies. Before this 
combined signal gets to the detector it still has to go through the Galactic magnetic field, where it 
is rotated once more by the amount RM^. The final measurement, RM^ is thus the linear sum 

RMm{z^) = RM} + mic{zi) + RM/(i;) + miG{zi) , (21) 

where zi is the redshift- vector [actually translated to a distance vector (c/. §2)] for the z*^ source. 
Since we are interested in the cosmological contribution to the RM [i.e. the second term in Eq. (^l])], 
the first step should involve "cleaning" the measured signal of all irrelevant extra contributions. 
We shall attempt to perform this cleaning in a statistical way. 
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3.2. 



The Galactic Mask 



There are two alternatives to assess the Galactic contribution to the measured RM. The two 
ways differ by the population sample used for the assessment. If an independent population at 
the outskirts of the Galaxy exists, for which RM can be measured, this population can be used to 
map out the Galactic RM. We shall hereafter use the term "mask source population" for the set 
of objects by which we map the Galactic contribution to the RM. 

An estimate for the thickness of the Galactic magnetic layer is ~ 1 kpc (Simard-Normandin 
& Kronberg 1980). That means that apart from the Galactic center direction, sources of distances 
that satisfy r^^c > 1/ sin(|6|) kpc (~ 3 kpc for |6| = 20°) are located beyond the Reynolds layer and 
fully probe the Galactic contribution to the RM. We hence consider two possibilities for the mask 
source population. 

One natural candidate for this role is the pulsar population, for which RMs are available. In 
order to make use of the pulsars we need to find a subset of them that reside in the outer part of 
the Milky- way magnetic layer. 

A class of pulsars that is especially appropriate for the task of probing the Galactic RM contribution 
is the millisecond pulsars. This population of old pulsars is believed to reside far out of the Galactic 
plane. The number density of millisecond pulsars as implied by a number of surveys at high Galactic 
latitude (~ 1 mJy sensitivity at ~ 1 GHz) ranges between 0.01-0.0175 deg-^, i.e. 400-700 pulsars 
over the sky (Foster, Cadwell, Wolszczan, &: Anderson 1995; Camilo, Nice, & Taylor 1996). 

The alternative mask population for assessment of the Galactic RM contribution is a sub-set 
of the closest extragalactic sources. The disadvantage of using this population is the need for a 
large enough source number at low enough redshift, to allow differentiation between the cosmological 
contribution and the Galactic one. Exploiting extragalactic sources for the resolution of the Galactic 
contribution affects only slightly the amount of cosmological contribution correlation, as long as 
the nearby sources are taken within z ^ 0.1, since the bulk of the sources is at z 1 (c/. §^). 

In order to use either one of the galactic mask populations, its number density should exceed 
a certain minimal number density. We now attempt to assess this minimal number. 
The Galactic RM contribution is a two-dimensional field with RM values. The data on the other 
hand are made of distinctive sources located at discrete directions. The lines-of-sight to the Galac- 
tic mask population, do not necessarily coincide with the line-of-sight direction towards the extra- 
galactic sources. The way to circumvent this difficulty is by introducing a smoothed version of the 
Galactic RM map, namely the field RMg'(z). The two fields are connected via the smoothing angle 
9s, and by a specific choice of a Gaussian smoothing for Nq measurements of the Galactic RMG(zj) 
located at the directions Zi 




i=l 



(22) 
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Fig. 3. — The 1° smoothed RM Galactic mask projected onto celestial coordinates. Positive RM values 
marked by continuous lines and negative values by dashed lines. The RM = line is thicker. The spacing 
is ARM = 30radm~^. The locations of sources in the SKB catalog (projected) are marked by solid dots, 
and the Galactic plane is marked by the dotted strip. 



When implemented using real data, the smoothing involves weights by measurement errors as well. 
For the purpose of this paper, we ignore these weights to avoid unnecessary complication. Using a 
smoothed field for the Galactic RM contribution, introduces random noise, ec, that can be directly 
evaluated by the very same Nq sample via 

'G = J^J2 (RMg(z,) - RMciz,)) . (23) 

The noise can be a function of position, it is difficult however to account for a position- 
dependent noise with low number of Nq- There are a few estimates in the literature for the noise 
level, ec, on different smoothing scales. 

An upper limit for ec is given in Minter & Spangler (1996), who calculate the observed "structure 
function", Drm which is the same as Eq. (^) but with an arbitrary lag 7 and the unsmoothed 
values of the RM. They conclude that L'rm(7 > 1°) = (340 ± 30)7°-''^=^°-°^rad2m-4, and drops 
steeper for 7 < 1°. The value of the Drm is an upper limit for €^(7) as calculated here, due to the 
smoothing we apply (and was not applied in the I^rm calculation). Simonetti, Cordes, & Spangler 
(1984) calculate Z?rm = 484 — 8100 rad^m"^, depending on the Galactic latitude, for much larger 
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Fig. 4. — The RM residuals from tlie 1° smootlied Galactic mask as function of the Galactic longitude (left) 
and latitude (right). 
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Fig. 5. — The 1° smoothed RM residuals of the Galactic mask in celestial coordinates. Symbols as in figure 



angular scales of 30° — 50° ! (linear scale). 

Simonetti & Cordes (1986) corroborate these results and consider even larger angular scales to 
obtain typically Urm < lOOOrad^m"^ on scales less than 100°. 
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Oren & Wolfe (1995) calculate a very similar quantity to ec (with a varying top hat window instead 
of a fixed Gaussian), and obtain a typical variance of < 1000 rad^m-^ on ~ 30° scales. 
All of these measurements suggest that for a mask population with number density of about one 
source per lOOOdeg.2, the Galactic contribution to the RM (|6| > 20°) can be resolved with a la 
accuracy of 30radm~^. That means that an isotropic coverage of about 100 — 200 mask sources 
on the sky is enough for achieving this noise level. 

"Damage control" for ignoring the position dependence can be devised by comparing the "clean" 
RM field dipole and quadrupole to the Galactic direction of the two. These first two multipole 
should be the most prominent signature of the Galactic RM mask. 

We go back to the suggested mask populations and examine whether they are able to fulfill 
the requirement of the minimal number density. 

The largest pulsar catalog to date [the electronic version of Taylor, Manchester, & Lyne (1993)] 
consists of 800 pulsars out of which 259 have RM measurements. This catalog was not compiled 
with a specific emphasis for millisecond pulsars, and therefore only 124 pulsars with RM values 
lie in Galactic latitude \b\ > 5°. The softer condition Vgrc > 2/sin(|6|) {i.e. considering an "RM 
layer" of 4 kpc) is fulfilled for 237 pulsars with both RM registered value and estimated distance 
r^rc (Taylor &z Cordes 1993). 

However, these 237 pulsars are still not distributed isotropically across the sky, and tend to con- 
centrate near the Galactic plane. 

A way to avoid this anisotropy is to compile a sample of millisecond pulsars, that lie far away 
from the Galactic plane. This population has a number density that exceeds the necessary ~ 10^^ 
deg.~^ by at least a factor of two and it may be as high as seven times denser than the minimal 
value. This population can be detected for RM (Fernando Camilo, private communication, 1997) 
and ensures the ec — 30 rad m~^ noise level. This kind of sample is ideal for the current paper 
analysis, as it does not contain any cosmological contribution or internal RM contribution. The 
main disadvantage of it is that it does not exist yet. 

In considering the extragalactic mask population; for the ^ 400 sources needed to resolve the 
Galactic contribution, it is sufficient to have a sample of ImJy sensitivity (at a few GHz) up to 
z ~ 0.1 [see Loan, Wall, & Lahav (1997) for the number density of radio sources with flux > 35 
mJy, and Dunlop & Peacock (1990) for the combined flat /steep spectrum luminosity function for 
extrapolation to the ImJy limit]. 

For this paper we follow Oren & Wolfe (1995) and use an extragalactic sample of RM measure- 
ments as an upper limit for the noise that appears by the Galactic RM removal procedure. The 
sample we use consists of 555 sources as listed by SKB. The SKB catalog allows us to realistically 
mimic the Galactic mask, and is an overestimate since we ascribe all the RM of the sources to the 
Galactic contribution. Moreover, unlike Oren & Wolfe we do not exclude "outliers" . This use aims 



^In the Taylor, Manchester, & Lyne catalog there are two pulsars with both RM registered value and p < 10 ms, 
and five for which p < 100 ms. The smallest polarized flux at 1.4 GHz for a pulsar with a registered RM is 0.2 mJy 
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to provide a realistic model of the Galactic mask contribution, and allows us to demonstrate the 
method power. For our purpose it should not be taken literally as the true Galactic mask contribu- 
tion (even though Oren & Wolfe did consider it that way) because no segregation by redshift was 
applied to choose mask sources from the SKB list. We neither attempted to evaluate the internal 
RM (by means explained in § p.3[ ) in order to minimize the ec value. 

The SKB source locations (projected, and therefore concentrated toward the circumference) 
are marked as black dots on Figure ^ and shown on top of the smoothed Galactic RM field with 
smoothing scale of 6s = 1°. The smooth map is shown in celestial coordinates. Figure ^ shows the 
RM residuals as functions of Galactic latitude and longitude, and figure ^ shows the contour map 
of the smoothed residual field. Both figures (|^ & ^) show hardly any longitude dependence of the 
residuals (in spite of the solar system off-center position), and bigger residuals in the ~ ±20° strip 
about the Galactic plane (c/. the same conclusion of Oren & Wolfe 1995). In an all-sky coverage 
of sources, this residual deviation should be of some worry, but as we shall see (p.2|), it doesn't 
bias the result for the most likely power spectrum as calculated by the Bayesian analysis. If the 
analysis is confined to high latitudes, then no systematic error due to the Galactic mask is expected 
(c/. For a few square degrees area, we notice there are typically no significant gradients in 
the Galactic RM field beyond \b\ ^ 30°. This allows us to use the smooth value safely when 
analyzing a small area with extragalactic RM measurements. This value is also consistent with all 
the abovementioned values of -Drm as obtained by more detailed analyses. 



3.3. The Internal Variation and foreground screens 

We would like to get an estimate for the internal RM contribution and for the foreground 
screen contribution (if the latter exists) . We do not attempt to correct each individual source for 
the RM contribution due to the internal and foreground screen terms. We do not take the measured 
RM and subtract an estimate for RM/+j, we rather attempt to estimate the distribution of RM/^j, 
and treat it as another noise term. 

The difference between the two separate terms in Eq. ( [2l| ) (RM/ and RMj ) emerges from 
the fact that while the internal contribution is obtained from the same medium that emits the 
polarized light, the foreground screen only serves as a filter for an already existing signal. In 
previous investigations, attempts were made to separate the two contributions (Burn 1966; Laing 
1984), later on it became clear that the signature of the two is rather similar (Tribble 1991). 

For our purposes we are interested in all previous calculations' "left overs", their signal is 
our noise and vice versa. For the noise estimate we claim it is legitimate to take the internal 
contribution and the foreground screen contribution together. We argue that the two can be dealt 
with simultaneously. To this end we shall use the other piece of information provided by the 
observations - the polarization degree. 

Tribble (1991) shows the connection between the observed depolarization (as a function of 
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wavelength) and the observed RM for an extended source. Using this connection, and provided 

the polarization measurements are available, a direct estimate for the internal RM can be carried 
out. Recall we do not expect any depolarization due to the cosmological RM since the cosmological 
coherence we are after is orders of magnitude bigger ('<^ a few h^^Mpc) than the source size. 

It is interesting to note, that even if we had perfect information about the depolarization 
due to the source structure and foreground screen contribution and an exact relation between the 
depolarization degree and the RM, we could still not subtract this derived RM from the observed 
one. This is due to the fact that the polarization degree doesn't specify the RM direction. 

Tribble's statistical connection between the observed depolarization and the observed RM is 
valid only if the correlation scale of the RM structure function for the source or the foreground 
screen is much shorter than the telescope beam size. This condition is probably not valid for 
individual damped Lyman-a systems, and galaxies that serve as foreground screens, and for which 
long range correlation across the telescope beam may exist. 

Full justification for neglecting damped Lyman-ct systems and galaxies exists only if we can 
avoid lines-of-sight that cross such systems. Otherwise the order of magnitude of an intervening 
galaxy contribution is obtained as follows: For a galaxy seen face-on, where there is a danger of a 

well ordered magnetic field, the disc thickness, to which the magnetic field is presumably connected 
is of the order of kpc (Simard-Normandin & Kronberg 1980). Observations show magnetic field 
magnitude of ^ /iG, mainly in the disc plane. In order to equate this to a cosmological magnetic 
field of the order of nG, and coherent scale of 50 h^^Mpc, the average baryon number density across 
the galaxy should be at least fifty times higher than the cosmological n^. 

However, we do not expect any correlation between the magnetic field orientation of galaxies 
along the line-of-sight, or galaxies adjacent (in angle) to each other. The worst contribution to 
the RM, could only come from a random walk of Nj steps where Nf is the number of foreground 
screens (intervening systems) along the line-of-sight. Since one can model the probability for an 
intersecting system along the line-of-sight as a function of the source redshift (Welter, Perry, &; 
Kronberg 1984; Lanzetta, Wolfe, & Turnshek 1995), another noise term can be added. We did not 
attempt to model this noise term here, because we believe that the best strategy would be to avoid 
highly contributing (i.e. damped Lyman-o; and Lyman limit) intervening systems. Low column 
density Lyman-ct system contribute very little to the RM signal (see Eq. (1.2) in Welter, Perry, & 
Kronberg 1984). An exception may exist, if two sources (or two separated parts of the same source) 
cross the same cluster. In that case the noise of the two sources (RM/_|_y) may be correlated, but 
then we can go back to the technique that uses the depolarization, in an attempt to remove the 
noise correlation. 

The best strategy, as stated earlier, would be to avoid intervening systems altogether. There 
are two methods for doing this. The first method is simply to avoid intervening systems by looking 
for sources that exhibit no absorption lines in their spectra (we assume spectrography is carried 
out anyway, to allow redshift determination). The existence of a large fraction of qiiasar lincs-of- 
sight which do not intersect a dense intervenor is corroborated by M0ller &; Jakobsen (1990) and 
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Lanzetta, Wolfe, & Turnshek (1995) analyses. 

The other method is to put a depolarization limit on the observations {e.g. the one imposed 
by Tabara & Inoue 1979). Such a limit naturally reduces the noise term of £/+/ and does not affect 
the correlation signal. As a matter of fact, a sequence of such depolarization limits, may be helpful 
in constituting the estimate for e/+/. There are yet more methods for minimizing the number of 
sources that manifest internal RM as we shall point out when we consider real observations (§^). 



3.4. Bayesian Likelihood Analysis 

In the current stage of research we have very little knowledge about the power spectrum of the 
cosmological magnetic field. We can hardly even put limits on its integral value, i.e. the rms value 
of the field. There are, however predictions regarding the power spectrum shape and amplitude. In 
the lack of any observational preference towards any of these predictions (apart from the limits on 
its rms value), we assume that all models are equally probable. Conventional estimates of the RM 
correlation, estimates of the sort applied to galaxies are not very useful. It is difficult to get correct 
error estimate without simulations, that in turn must assume some magnetic field power spectrum, 
moreover in the traditional correlation calculation procedure, the errors of different bins in r space 
are correlated. On top of that, the quantity we are after is the magnetic field power spectrum, 
and the inversion from the integrated RM back to the magnetic field is a non-trivial one. All of 
the above lead us to the use of Bayesian statistics as a tool for finding the best parameters, and 
their probability, given a model. In the Bayesian formalism, the a-posteriori probability density of 
a model, m, given the data, d, is 

^Hd) = «^. (24) 

As stated earlier, P(m), the model probability density is unknown and therefore assumed uniform 
for all models. The data probability, in the denominator is the same for all models, and therefore 
serves as a normalization factor. We are thus left with equivalence between the operation of 
maximizing the probability of the model, given the data [P(m|d)] and the operation of maximizing 
the probability of the data, given the model [P(d|m)]. This equivalence allows us to write down 
the likelihood function for Ngrc sources with measured RM. 

We begin by constructing the Ngrc x ]^src symmetric matrix Tij of the expectation values for 
the RM correlation (Eq. [l^) between these sources (as determined by their position). We further 
assume that the errors in each RM value are uncorrelated and Gaussianly distributed. The variance 
of the overall error is the sum of the quadratures of the various error sources i.e. 



e: 



and the likelihood function becomes 

N'src 



C = [(2^)^- det(Ti,)]-V2 exp f -1 ^ RM.Tr.iRM,- | , (26) 
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where fij = Tij + 

The statistic is defined as = — 21ni2. The statistics as defined, is a distribution (of 
Ngrc degrees of freedom) with respect to the data points, not the model parameters. This x^ is 
interpreted as a rough estimate for the confidence levels in the parameter space. 



3.5. Estimating the Signal-to-Noise Ratio 

Before proceeding to the elaborated description of the various tests we performed with artificial 
data, we would like to get a rough estimate for the expected RM correlation and the necessary 
number of sources for a high enough signal-to-noise ratio at various scales. We temporarily abandon 
the Bayesian approach, and focus on V{m\d). Recall that the noise in the cross correlation is not 
correlated, but may still dominate the signal if not appropriately averaged over many pairs. The 
signal-to- noise ratio for an individual measurement of a pair at identical redshift, a separation angle 
7, and in the context of a specific model, is given by 

T(z,7)\_ (T(z,7)) ^27) 



£2 / (e2) ' 

where in the last equality we assume no dependence between the various noise terms and the model 
(cosmology -|- power spectrum). The cosmological parameters we choose to use are as follows: We 
restrict ourselves to fiat cosmology with 0.^ = Ij we choose Hq = 100 kms~^Mpc~^ (but the 
scaling with h is straightforward). We take the ionization factor to be Xn = 1 which according 
to the Gunn-Peterson effect (Gunn & Peterson 1965) is a reasonable choice, and = 0.024/1"^ 
(Tytler, Fan, & Buries 1996), i.e. fib = 2.6 x lO"'^ cm'^ 

For Npair pairs where each source is a member of one pair only, the signal-to-noise can further 
be increased by y/ Npair ■ In figure |^ we plot the necessary number of pairs in order to achieve a 
signal to noise ratio of Q = 3. All sources are assumed to be at the same redshift (2, 1, or 0.5), 
we take the noise values ec = 32radm~^ (this is the value for 6s = 1°). For the noise term due to 
the internal RM and the (possible) foreground screen we use the Tabara & Inoue (1979) catalog as 
a guideline. Tabara Sz Inoue listed ~ 1500 radio sources in their catalog with RM ^ 200radm~^. 
That means a ~ 3.4a cover, and gives (JRM,obs — 59radm~^. The most conservative assumption 
one can make is to attribute all of this RM to the internal and foreground screen contribution. 
The procedure for deriving the internal and foreground screen contribution from the depolarization 
(c/. § ^) bears noise itself and is correct only statistically. We model its accuracy as a Gaussian 
with one quarter the width of the observed RM distribution. Doing so we end up with e/^j = 
(59^ -|- 15'^)^/2 ~ 61radm~2. This is a conservative choice, and the true value for £/+/ is probably 
much smaller. Note for instance the £/+/ — 20 radm"^ value as obtained by Oren &: Wolfe (1995). 
The measurement errors are neglected {of. § 3.1). Two power spectra are considered and they are 



both normalized to give BJ'J^^{50 h-^Mpc) = 1 nG (the n = case gives results between the n = — 1 
and the n = 1 power spectra). A typical number of pairs for achieving a signal to noise ratio of 
~ 3 in the range 7 ^ 0.5 is between 10 and 2 x 10'^ pairs, depending on the power spectrum. 
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Fig. 6. — Estimate of the number of pairs at separation 7 needed to achieve signal to noise ratio of 3. The 
magnetic field is normalized by .6^^(50 h~^Mpc) = 1 nG, and the noise calculation is explained in the 
text. Two cases for power index of —1 (left) and 1 (right) are plotted. 

This estimate of Q is not accurate. To begin with, for Ngrc sources we expect less than 
{Ng^^ — Ngrc)/'^ statistically independent pairs. We do not expect all sources to reside at the same 
redshift, neither do we expect the errors to be identical for all sources. Furthermore, the separation 
between the sources is not fixed , but rather is spread from the arc-minute scale to 180°. Due 
to the need for a realistic error estimate, and a demonstration of an actual implementation of the 
technique we propose, led us to apply it to mock catalogs, for which the underlying power spectrum 
is known. 



4. TESTING WITH ARTIFICIAL DATA 

As a reliable probe for our method we create artificial catalogs of radio sources for which RM 
values are measured. We first generate a Gaussian random field of the magnetic vector potential, 
A with the power spectrum PA{k) = k'^PB{k){= E{k)). The realization takes place on a 256^ grid 
points equally spaced in comoving coordinates. The Psik) power spectrum is Gaussianly smoothed 
on the one grid cell scale. The amplitude of the field is arbitrary. We then derive the magnetic field 
itself by the real space relation ^ = V x ^ translated to k space i.e. B{k) = —ik x A. We use the 
k space symmetries due to the real nature of the A field (no imaginary part), and go back to real 
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space to obtain a divergence-free magnetic field with the desired correlation function, and periodic 
boundary conditions. We check the resultant field by calculating the field divergence about each 
and every grid point in real space and get V • B/B^ms < 5 x 10^^. 

At this point we select the cosmology into which we embed the simulation. Our standard 
choice is identical to the one described in the signal-to-noise estimate (§ |3.5| ). Due to the limited 
dynamical range, we choose the grid scale to represent 1.5 h~^Mpc, and we are therefore confined 
to a largest simulation wavelength of 768 h^^Mpc. 

We then select the source mask, i.e. the source distribution across the sky, and their redshifts. 
The two-dimensional distribution is taken to be a Poisson-like distribution, this choice allows us 
to have close pairs (unlike grid selection), but to be as conservative as possible in terms of the 
angular correlation function of the sources. Introducing any correlation as observed for radio 
sources (Kooiman, Burne &: Klypin 1995; Sicotte 1995; Cress et al. 1995; Loan, Wall, &: Lahav 
1997) can only increase the number of close pairs on small angles, and thus improve the RM 
correlation resultant signal (see fig. ^). We do not take explicitly into account very small angle 
(sub arc-minute) separations that may exist in real data sets for resolved extended sources, whose 
different parts can be used as more than one background source. We employ two distribution 
schemes, an all-sky coverage, and a cover of 150 deg.^ area centered about the north Galactic pole. 

The source redshift distribution is selected according to the N{z) of radio sources with flux 
> 35 mJy (at 4.85 GHz) as calculated by Loan, Wall, &: Lahav (1997) who used the mean of 
the theoretical luminosity function models of Dunlop & Peacock (1990) and kindly provided their 
fit. Loan et al. (1997) N{z) is peaked around z = 1, and has a shape resembling a low redshift 
truncated Gaussian of width (2 x a) Az ~ 1.9. 

In principle, setting flux limit and sky coverage, sets the source number. In practice only a 
sub-set of all radio sources emit polarized light. A random sample taken from the NVSS survey 
(Condon et al. 1994) shows registered polarization angle for ~ 40% of the sources. Refraining from 
sources with intervening systems and large internal RM may further decrease the number. Since 
the source number density in the radio catalogs of this fiux limit is a few per square degree, when 
we prepare an "all-sky" coverage catalog (see below), we are always way below the available number 
of sources. For smaller angular coverage with higher source density, one has to further decrease 
the flux limit. For example in the 2.5 mJy limit we expect ~ 50 sources per square degree (NVSS, 
Condon et al. 1994), and in the 1 mJy limit we expect source density of about 100 per square 
degree (the "FIRST" survey; Becker, White, & Helfand 1995). 

The N{z) function of Loan et al. (1997) is therefore interpreted as a normalized selection 
function, with a cutoff at the highest redshift of the catalog. We expect this selection function to 
provide an underestimate for the number of sources at high redshift, if a lower fiux limit is set. 
Since sources at higher redshift typically bear higher signal (without changing the noise), such an 
imposed selection is a conservative choice in terms of the expected signal-to-noise ratio 0. 



^The N{z) we choose, which is peaked in a relatively low redshift, may be useful in mimicking another effect. As 
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For each selected source, we integrate along the line-of-sight over the dot product between 
the magnetic field and the line-of-sight direction using the periodic boundary conditions. The 
integration scheme assigns weights to each integration segment, k, according to 

^^-Krlizk) 
in order to account for the cosmological evolution. 

We then add the following noise terms to the resultant RM 

1. Internal variation + foreground screens. 

We draw this noise term from a Gaussian of (59radm^^)^ variance (c/. § 3^). This is the 



added RM to the cosmological one. However, we also try to imitate the procedure of e/+/ 
evaluation by the depolarization. To this end, for a specific value of RM}_,_j (drawn earlier), 
we mimic the recovery of the £/+/ from the depolarization degree, by listing an £/+/ value 
scattered about RM}^^. The scatter has a Gaussian distribution of (15radm~^)^ variance. 
This is the error in the noise estimate. The "observer" takes into account the scattered value 
of the and not the value that was actually added to the integrated cosmological RM. 

2. Galactic mask. 

We use the model-estimate of the Galactic mask, smoothed to 1° scale (see § |3.2| ). Since the 
smoothed galactic RM differs from the true RM to the line-of-sight, we add a RM drawn from 
a Gaussian distribution with the standard deviation of ec = 32radm^^ (Eq. p^). 

At the end of this process we are left with a mock RM catalog that consists of source coordinates 
and redshift, measured RM value, and the total error in the RM value {i.e. ^i+f). 

We proceed by feeding the mock catalog into the maximum likelihood procedure. In the 
current application, we consider only the right cosmology, i.e. the one in which the simulation was 
embedded. In principle the sensitivity to the cosmology choice can be checked as well. The two 
variables for the model are the amplitude and the power index. It is really the combination of the 
two that the likelihood procedure constraints best. 

The statistical significance of the distribution (as defined following Eq. (^)) is obtained in the 
standard way. The number of degrees of freedom is the source number minus the number of fitting 
parameters, and the goodness of fit is calculated by taking into account the value of Xmin- 
goodness of fit turns out to be very small (< 0.1), this is a clear indication for convergence to the 
wrong minimum, and vice versa. For the right solution the minimum value for the P^r degree 
of freedom (the source number) is typically unity within 10~^ accuracy. 



the source redshift increases, the probability for intervening systems along the line-of-sight increases as well. If in a 
sample compilation we try to avoid intervening system, this attempt translates to a steeper fall-off of the selection 
function. 
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4.1. Partial Sky Coverage 

We first explore partial sky coverage mock RM catalogs. Figures ^ show contour maps of 
— Xmin with steps of unity. The thicker lines indicate the approximate 1,2, and 3a levels of 
the probability function as ascribed by using two parameter statistics (c/. § |3.4|) . The power 
spectrum normalization is our standard BjJ^g{50 h~^Mpc) = 1 — 5 nG, depending on the power 
spectrum. We use 500 — 800 sources for the part sky and 150 square degrees coverage. The different 
source number and Brms values reflect the fact that we have fixed the angular sky coverage. That 
in turn translates differently for different spectra. A redshift cutoff of z < 5 has been applied. 
Figure § shows a particular example for the case n = 1, where we have increased Bj^g{50 h^^Mpc) 
to the 17 nG level in order to verify the result dependence on Brms- 

We notice that for all power spectra, the Bayesian approach recovers the correct value from 
the artificial catalog (marked by a cross sign) to within the Icr level. The method is more robust 
with respect to the amplitude determination. The power index is not well bound unless it is very 
different from the true one. The method (with the limited number of sources) tends to recover best 
the Brms value on a certain scale, which is a combination of the power index and the amplitude. 
Typical values for the ratio between the true Brms and the minimization result for Brms are in the 
range 0.6 — 1.5. 

Partial sky coverage is preferable in terms of the number of small separation pairs, and the 
indifference to the Galactic mask variation as function of direction. A collection of several of 
these small patches is expected to provide an even better constrain on the magnetic field power 
spectrum, but the analysis becomes a bit more complicated, as one should also include the large 
scale correlation between the patches without increasing the noise from the auto correlation to an 
intolerable level. 

Given the relatively big errors, we desire off-diagonal terms in the correlation matrix to be 
as big as possible {i.e. small separation). This is meant to avoid heavy weight for diagonal terms 
which are dominated by noise. The price for choosing small patches with small sky coverage is 
limited ability to constrain the power index. Under special conditions, this ability can be partially 
recovered by turning to an all sky coverage as we discuss in the next section. 

4.2. All Sky Coverage 

An all sky coverage, with the small correlation scales that are tested here, is impractical without 
lowering the error level significantly. Otherwise the Bayesian analysis will all be dominated by the 
diagonal terms with a large noise added to them and one loses the advantage of the correlation 
scheme for noise reduction. 

All sky coverage is however still feasible if we design the sample very carefully. In order to 
minimize the Galactic noise we should observe sources in proximity to the mask sources direction 
and avoid the necessity for smoothing the Galactic mask. Instead, the RM values as deduced 
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Fig. 7. — Contour plot of ~Xmin (^^ defined following Eq. (p6|)) in the amplitude - spectral index plane. 
The amplitude is in units of the true simulation amplitude. The spacing is of = 1, and confidence 

levels for 1,2, and 3a are plotted (bold lines) as if the contours obey statistics in the two-dimensional 
parameter space. The true value of the simulation is marked by the cross sign. All plots are for a part sky 
coverage of 150 deg.^ about the Galactic north pole. The values for BjJ^g{50 h~^Mpc) are 5, 2, and 1 for 
n = —1, 0, and 1 (top, left, right) respectively. The total number of sources is marked on the plots. 

from the mask source population can simply be subtracted from the observed RM values. The 
small variation in the Galactic magnetic field on small scales (Simonetti, Cordes, & Spangler 1984; 
Simonetti & Cordes 1986; Minter & Spangler 1996) allows this subtraction to take place. 

If on top of that one can reduce the noise from internal RM to the level of 5 radm~^ (one 
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Fig. 8. — Same as the n = panel in the previous figure (fig. 0). Here Bj^g{50 h -"^Mpc) is 17 nG, and 
allows a better determination of the parameters. 



can always avoid foreground screens) by multi wavelength observations, then it becomes sensible 
to exploit an all-sky coverage. 

Figure |9| shows an example of a sample taken at the mask source directions with ej = 5 radm~^ 
and Bjj^g{50 h~^Mpc) = 1 nG. The rest of the parameters are identical to those of the part 
sky coverage (n = 1) of the last section. We notice that the power index is somewhat better 
constrained in an all sky coverage, especially if the amplitude (or Brms) are assumed or known 
from some different source. The true values are recovered to within the ~ 2a level. Moreover 
such coverage allows us to calculate all prevailing magnetic fields on the sample scale, expressed as 
Pb(/c) = S^{k — ^) with 6^ the Dirac delta function and R the sample depth. Notice, however 
that using the formalism proposed in section ^, we cannot take into account power spectra that are 
not functions of |A;| alone so this test is a bit weaker than the dipole test. 



5. DISCUSSION AND CONCLUSIONS 

We have presented the formalism for the RM correlation function and demonstrated that by 
using it we can put limits on the power spectrum of the cosmological magnetic field. These limits 
are more stringent than limits obtained by any other method that currently exists. Limits of 2 — 3cr 
level for magnetic field of ~ 1 nG on inter-cluster scales (~ 50 h~^Mpc) can be devised with only 
10^ — 10^ sources. The increment of the source number reflects linearly in the Brms limits. 

We performed the statistical analysis following a Bayesian approach that seems to be the most 
adequate for this problem. The correlation method we propose is also less susceptible to systematic 
effects such as observational bias (or bias due to evolution) toward smaller internal RM measures 
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Fig. 9. — Contour plot of ~ Xmin ™ amplitude - spectral index plane (same as figure ^) for an 
all-sky coverage of sources in the mask sources directions. The internal scatter in RM is assumed to be 
considerably lower than in the part sky coverage and 5^^,(50 h-^Mpc) = 1 nG. 

for higher redshift sources. Both of these biases tend to lower the expected correlation signal as 
the source redshifts increase, whereas a cosmological magnetic field tends to raise the correlation 
signal between such sources. 

In order to estimate the statistical significance of limits as derived by the proposed method, we 
have simulated a few possible cosmological magnetic fields. The simulations allow us to have a fair 
estimate of the errors involved in the power spectrum derivation. The analysis depends, however 
on a few assumptions, and most crucially on the ability to eliminate or estimate the foreground 
screen and internal contribution to the observed RM. 

This crucial point can be addressed in the process of the sample compilation. Reduction of 
the foreground screen and internal RM contribution is achievable by applying two cutoff criteria 
to the sample. The first cutoff in minimal polarization degree is meant to avoid sources that 
exhibit large depolarization degree. These sources are usually interpreted as having a big non- 
cosmological contribution to their measured RM. The second cutoff is for sources that do not follow 
a dependence of the RM. Following Laing's (1984) recipe, in order to minimize the number of 
sources with internal RM, the wavelength span should allow detection of deviation from the 
dependence of the RM for (p > 7r/2. Should this deviation occur, the source is to be eliminated 
from the catalog as a suspect for substantial internal RM contribution. 

On top of the abovementioned cutoffs, all sources should preferably be detected for absorption 
lines, as a method to eliminate sources with intervening systems along the line-of-sight to the 
sources. 

The spatial coverage of the sky should be twofold. Small patches of a few square degrees 
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enable the small separation pair number needed for noise reduction. Small patches, however, can 
not reveal large correlation scales in an effective way. It is thus desirable to compile a dilute all sky 
coverage sample, in addition to the high number density patches off the Galactic plane. The two 
sampling strategies cover a large range of potential correlation scales, without the necessity for an 
all-sky dense sample. The all sky coverage should preferably look for sources in proximity to the 
mask sources directions in order to minimize the Galactic noise term. Accurate internal RM for 
these sources should be evaluated through depolarization measurements in order for them to be 
useful for the power spectrum evaluation. Ultimately, the sample can always be mimicked in terms 
of the exact source locations in the framework of all and every detected model. Exact imitation of 
the source locations ensures the right correlation terms in the correlation matrix. 

The limits we have derived here are a combination of limits on the free electron average density, 
fie, and the magnetic field. If by any other fashion (like HII absorption) an estimate for can 
be achieved, then the limits on B will become more robust. The estimate is also a function of the 
assumed cosmology, and may be entangled with magnetic field evolution in the post recombination 
era. We haven't modeled such evolution in this paper. 

This method of RM correlation can further be generalized to the smoothed RM correlation. 
The data can be smoothed on a certain angular scale, and then either compared to the expectation 
value predicted by a model, or inverted numerically to give limits on B. The Bayesian approach, 
though ceases to be advantageous in the smoothed case, because even though the noise terms 
decrease, the smoothing introduces correlation among the errors, and the statistical interpretation 
becomes more hazardous. 

A nice feature of the proposed analysis is that it is bound to yield results. These can be either 
limits on the magnetic field magnitude, or actual detection of its value. Either way, these results 
may help in lifting the curtain over the mystery of the cluster and galactic magnetic fields origin. 
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